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GENERAL  DISCUSSION 
In  the  following  discussion  we  will  summarize  the 
work  performed  by  Science  Applications,  Inc.  (SAI)  on 
"Energy  Release  and  Fluid  Dynamics  in  Multiphase  Systems" 
(Contract  #N00014-82-C-2037 ,  SAI  Project  #1-157-18-519-00) 
during  the  technical  performance  period  ending  7  January 
1983.  During  this  contract  period,  our  attention  was 
focussed  primarily  on  (1)  deriving  power-energy  relations 
for  the  direct  initiation  of  detonations,  and  (2)  ignition 
behind  incident  shocks  and  the  transition  to  detonation. 

The  efforts  and  accomplishments  in  each  of  these  areas 
described  below  in  some  detail. 

Power-energy  Relations  for  the  Initiation  of  Detonations 

i  o  7  j 

Recent  experimental  studies  *  *  on  the  direct 

initiation  of  gaseous  detonations  have  shown  that  initiation 

depends  not  only  on  the  energy  deposited  but  also  on  the 

rate  at  which  it  is  deposited,  namely  the  power.  The  experi- 
2 

mental  results  have  also  indicated  that  there  is  a  minimum 
initiation  energy  below  which  a  detonation  will  not  occur 
no  matter  what  the  power  is  and  that  there  is  a  minimum 
power  below  which  a  detonation  will  not  occur  no  matter 
what  the  total  energy  is.  The  relation  between  the  power 
and  the  energy  obtained  from  different  experimental  arrange¬ 
ments  have,  however,  exhibited  different  qualitative 
characteristics . 
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We  have  modified  and  extended  a  basic  theoretical 
model  proposed  by  Abouseif  and  Toong^  and  have  used  it  to 
determine  the  relation  between  the  power  and  the  energy 
required  for  the  initiation  of  planar,  cylindrical  and 
spherical  detonations  in  a  gas  mixture.  The  model  considers 
the  generation  of  a  constant  velocity  shock  by  means  of 
appropriate  energy  addition.  The  flow  fields  behind  the 
shocks  in  the  three  geometries  are  obtained  by  solving  the 
one-dimensional  conservation  equations  for  mass,  momentum 
and  energy.  The  solution  procedure  is  simplified  by  using 
a  similarity  solution.  The  model  also  requires  knowledge 
of  the  duration  for  which  the  shock  must  be  maintained  in 
order  to  initiate  a  detonation.  This  time  is  related  but 
not  necessarily  equal  to  the  chemical  induction  time  of 
the  gas  mixture. 

The  results  from  the  model  show  that  the  qualitative 

4 

differences  in  the  power-energy  relations  of  Dabora  and 
Knystautas  and  Lee^  are  due  to  the  differences  in  the 
geometry  of  the  detonations  in  the  two  cases.  The  results 
also  show  that  the  minimum  power  and  the  minimum  energy  are 
very  sensitive  to  uncertainties  in  the  ignition  delay  time 
used.  Another  interesting  observation  from  the  model  is 
that  the  minimum  power  required  corresponds  to  a  shock  of 
minimum  Mach  number  only  in  the  case  of  planar  detonations. 
In  the  cylindrical,  and  spherical  cases,  it  is  possible  to 
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initiate  a  detonation  with  a  shock  wave  of  lower  Mach 
number  than  that  corresponding  to  the  minimum  power.  Such 
a  shock  wave  will  have  to  be  maintained  for  a  longer  time 
than  the  shock  corresponding  to  the  minimum  power  and  hence 
will  require  a  larger  amount  of  energy. 

This  work  will  be  presented  at  the  Spring  Meeting 
of  the  Western  States  Section  of  the  Combustion  Institute 
as  Paper  No.  WSCI-83-35.  A  copy  of  this  paper  is  included 
in  this  report  as  Appendix  A  entitled,  "Power-Energy  Re¬ 
lations  for  the  Direct  Initiation  of  Detonations  in  Oxy- 
Acetylene  Mixtures."  A  more  detailed  account  of  this  work 
is  currently  being  prepared  for  publication  as  a  NRL 
Memorandum  Report. 

An  important  extension  of  the  work  described  above 
is  to  the  study  of  condensed  phase  mixtures.  For  this 
purpose  we  will  have  to  replace  the  perfect  gas  equation 
of  state  which  is  used  in  the  above  model  with  appropriate 
equations  of  state.  Furthermore,  instead  of  the  induction 
time  we  need  to  use  the  run  time  to  detonation  which  has 
been  measured  for  several  explosives^.  With  these  modifi¬ 
cations,  which  are  currently  being  pursued,  this  theoretical 
model  will  be  a  valuable  tool  in  classifying  condensed  phase 
explosives  according  to  their  sensitivity  to  detonations. 


Detonation 


i  o  9  in  ii 

Experiments'*  ’  ’  and  numerical  simulations  have 
shown  that  shock  induced  ignition  in  gaseous  hydrogen-air 
mixtures  may  occur  in  one  of  two  distinct  modes,  depending 
on  the  thermodynamic  state  in  the  shocked  mixture.  At  lower 
temperatures  the  ignition  is  weak,  or  mild,  with  the  gradual 
development  of  the  gas  dynamic  explosion;  at  higher  tempera¬ 
tures  it  is  strong,  or  sharp,  with  a  practically  instantane¬ 
ous  appearance  of  a  secondary  shock  induced  by  the  explosive 
reaction.  Furthermore,  at  low  temperatures  the  formation  of 
distinct  flame  kernels  appears  to  be  an  essential  precursor 
of  ignition. 

We  have  performed  time -dependent  numerical  simulations 

to  elucidate  some  of  the  details  of  weak  ignition  behind 

incident  shock  waves  and  the  subsequent  transition  to 

detonation.  These  simulations  were  carried  out  using  the 

12  13  - .  ‘ 

NRL  ID  reactive  shock  code  *  which  comb.ines  a  description 
of  the  fluid  dynamics  and  detailed  chemical  kinetics. 

The  numerical  simulations  show  that  a  small  amount  of 
energy  released  in  the  shocked  region  can  be  the  origin  of 
pressure  waves  which  accelerate  the  shock  front.  This 
shock  acceleration  leads  to  an  increase  in  the  temperature 
and  pressure  behind  it.  In  the  weak  ignition  regime,  a 
small  increase  in  the  pressure  and  temperature  can  result 


in  a  substantial  decrease  in  the  induction  time.  The 
numerical  simulations  show  how  this  leads  to  the  formation 
of  reaction  centers.  The  formation  of  a  hot  spot  due  to 
energy  release  at  one  of  the  reaction  centers  and  the 
subsequent  development  of  a  pair  of  reaction  waves  from 
the  hot  spot  have  been  studied.  The  acceleration  of  these 
reaction  waves  into  detonation  waves  has  also  been  observed 
in  the  simulations. 

The  predictions  of  the  model  are  in  qualitative  agree¬ 
ment  with  experimental  observations.  However,  the  model 
does  not  include  multi-dimensional  phenomena  such  as  tur¬ 
bulence  and  boundary  layer  growth  which  play  an  important 
part  in  any  quantitative  study  of  the  transition  to  detona¬ 
tion.  Currently,  a  two-dimensional  reactive  shock  model 

exists  but  it  uses  a  parameterized  model  for  energy 
14  15 

release  ’  .  The  model  is  now  being  extended  to  include 

a  detailed  chemical  kinetic  scheme.  Calculations  such  as 
discussed  above  performed  with  this  model  will  show  the 
effects  of  transverse  waves  and  boundary  layers  on  the 
transition  to  detonation. 

This  work  has  been  presented  at  different  stages  of 
completion  both  at  the  35th  Meeting  of  the  American  Physi¬ 
cal  Society,  Division  of  Fluid  Dynamics  and  at  the  1982 
Fall  Meeting  of  the  Eastern  States  Section  of  the  Combustion 
Institute.  A  short  version  of  these  presentations  is  included 
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in  this  report  as  Appendix  B  entitled,  "Numerical  Simulations 
of  Weak  Ignition  Behind  Incident  Shock  Waves  and  the  Transi¬ 
tion  to  Detonation".  A  more  detailed  and  complete  version  of 
this  work  is  presented  in  Appendix  C  entitled,  "Ignition  of 
Flamelets  Behind  Incident  Shock  Waves  and  the  Transition  to 
Detonation'1  and  this  has  been  accepted  for  publication  in 
the  journal.  Combustion  Science  and  Technology. 
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ABSTRACT 


Recent  experimental  studies  on  the  direct  initiation  of  gaseous 
detonations  have  shown  that  initiation  depends  not  only  on  the  energy 
deposited  but  also  on  the  rate  at  which  it  is  deposited,  namely  the  power. 
In  this  paper,  we  have  used  a  theoretical  model  to  determine  the  relation 
between  the  power  and  the  energy  required  for  the  initiation  of  planar, 
cylindrical  and  spherical  detonations  in  a  detonable  gas  mixture.  The 
results  from  the  model  show  that  the  qualitative  differences  in  the  power- 
energy  relations  obtained  from  two  different  experimental  arrangements 
are  due  to  differences  In  the  geometry.  We  also  show  that  the  minimum 
power  requirement  corresponds  to  a  shock  of  minimum  Mach  number  only  in 
the  case  of  planar  detonations.  Finally,  the  effect  on  the  power- 
energy  relation  of  the  experimental  uncertainties  in  the  determination 
of  the  induction  times  has  been  studied  for  a  stoichiometric  oxy-acetylene 


mixture 
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Introdoctlon 

The  early  stadias  of  dirset  initiation  of  gaseous  detonations1 ' *' 5 
established  the  importance  of  the  magnitude  of  the  sooroe  energy.  More 
recent  experiments**'  5'  8  here  shown  the  isportanoe  not  only  of  the  energy 
but  also  of  the  rate  at  whieh  the  energy  is  deposited,  namely  the  power.  3he 
experimental  results  of  Lee  et  al.  5  indicate  that  there  is  a  minium  detona¬ 
tion  energy,  Eg,  below  which  a  detonation  would  not  occur  no  matter  what 
the  power  is  and  that  there  is  a  minimum  power,  Eg,  below  which  a  detona¬ 
tion  would  not  occur  no  matter  what  the  total  energy  is.  later,  they  noted  6 
that  the  requirement  for  a  minimum  value  for  the  power  of  the  souroe  indi¬ 
cates  that  the  source  amst  be  capable  of  generating  a  shock  wave  of  certain 
minimum  strength  (Mach  nmaber).  They  also  concluded  that  the  minimum  energy 
requirement  isplied  that  the  shock  ware  mist  be  maintained  at  or  above  this 
strength  for  a  certain  minimum  duration. 

Aecently  these  ideas  have  been  used  by  Debora7' ®  to  obtain  a  relation 
between  the  power  and  energy  required  for  the  direct  initiation  of  hydrogen- 
air  detonations  in  a  shock  tube.  However,  this  power-energy  relation  is  very 


different  qualitatively  from  those  of  Xnystautas  and  Lee6.  More  recently 
Abouseif  and  Soon?*  have  proposed  a  sixple  theoretical  nodal  to  determine  the 
power-energy  relation  and  predict  their  respective  threshold  values.  3he 
predictions  based  on  their  model  were  in  qualitative  agreement  with  the 
experiments  of  Knystautas  and  Lee6. 

In  this  paper  we  have  modified  and  extended  the  basic  model  proposed  by 
Abouseif  and  Toong9  and  have  used  it  to  determine  the  relation  between  the 
power  and  the  energy  required  for  the  initiation  of  planar,  cylindrical  and 
spherical  detonations  in  a  datonable  gas  mixture.  Specifically,  we  discuss 
its  application  to  a  stoichiometric  oxy-aeetylene  mixture. 

The  Theoretical  Model 

The  model  considers  "the  flow  generated  by  the  motion  of  a  constant 
velocity  shock  wave  in  planar,  cylindrical  and  spherical  geometries.  As  this 
shock  wave  passes  through  a  gas  mixture,  the  gas  temperature  and  pressure 
increases.  Cue  to  this  increase  in  temperature  and  pressure,  ignition  could 
occur  in  the  shock  heated  gas  mixture  after  the  elapse  of  a  certain  time  and 
this  may  lead  to  a  detonation. 

A  constant  velocity  shockwave  can  be  formed  in  each  of  the  three 
geometries  by  the  motion  of  a,  constant  velocity  piston10'11.  Furthermore, 
it  has  been  shown12  that  a  pressure  and  velocity  field  identical  to  that 
ahead  of  a  constant  velocity  piston  can  be  generated  by  appropriate  heat 
addition.  Tor  exa^le,  a  flow  field  bounded  by  a  constant  velocity  planar 
piston  and  a  constant  velocity  planar  shock  wave  can  be  generated  by  a  planar 
energy  source  with  a  constant  rate  of  energy  deposition.  An  exasple  of  such 
an  energy  source  is  the  high  pressure  driver  in  a  constant  area  shock  tube. 
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In  general/  the  source  power  P  to  generate  a  constant  velocity  piston  In 
planar/  cylindrical  and  spherical  geometries  can  be  written  as9'*1 

“  T^iT  Vpup°  ta_1  '  (1> 

where  CQ  -  1,  2%,  4s  for  a  m  1/2,3  corresponding  to  the  planar,  cylindrical 
and  spherical  geometries  respectively;  pp  and  Up  are  the  pressure  and 
velocity  at  the  piston  surface  and  t  is  the  duration  of  energy  deposition* 
The  energy  deposited  is  given  by  the  time  integral  of  the  power,  that  is 


E(t) 


MM)  8 


(2) 


In  the  planar  case,  the  pressure  and  fluid  velocity  at  the  piston 
surface  are  the  same  as  those  just  behind  the  shock.  However,  in  the 
cylindrical  and  spherical  cases,  the  flew  field  between  the  shock  and  the 
piston  surface  is  nonuniform  and  can  be  obtained  by  solving  the  governing 
partial  differential  equations*  Bowever,  the  solution  procedure  is 
considerably  simplified  if  we  seek  a  s  ini  lari  ty  solution*  Then  the  system  of 
partial  differential  equations  can  be  reduced  to  a  system  of  coupled  ordinary 
differential  equations: 


dP  ,  fp  dp 
dX>  p  dL 


(5) 
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Zb  the  above  system  of  equations/  the  density  p,  the  velocity  u  and  the 

pressure  p  are  all  functions  of  the  similarity  variable  L,  which  is  equal  to 

the  radial  location  r  divided  by  the  time  t.  The  pressure  and  the  velocity 

at  the  piston  surface  which  are  required  in  Eqs.  (1,2)  can  be  obtained  by 

solving  Bqs.  (3-5)  in  the  following  manner,  for  a  shock  of  a  given  Mach 

number,  we  can  calculate  the  flow  condition  just  behind  the  shock  using 

normal  shock  relations.  We  can  then  integrate  Bqs.  (3-5)  from  just  behind 

the  shock  to  the  oiston  surface  to  obtain  p„  and  u  which  are  needed  in 

P  P 

Bqs.  (1,2).  The  procedure  is  further  sisplified  by  appropriately  combining 
Eqs.  (3-5)  into  two  equations  and  normalizing  them.  This  is  discussed  in 
detail  elsewhere11. 

In  order  to  determine  the  power-energy  relation  using  Bqs.  (1,2)  we  also 
need  to  know  the  duration  for  which  the  energy  mist  be  deposited  in  order  to 
initiate  a  detonation.  This  tine  must  at  least  be  equal  to  the  time  at  which 
ignition  first  occurs  in  the  flow  field9.  As  noted  by  Qrtiew  and  Cppenheim13 
ignition  usually  occurs  first  at  the  contact  surface  (i.e.,  at  the  piston 
surface  here)'  since  the  tesperature  and  pressure  is  highest  at  this  location. 
So  a  first  estimate  of  the  time  t  in  Bqs.  (1,2)  would  be  the  induction  delay 
time  corresponding  to  the  conditions  at  the  piston  surface. 

Results  and  Discussion 

We  first  used  the  model  described  above  to  determine  the  power-energy 
relation  for  the  initiation  of  cylindrical  detonations  in  a  stoichiometric 
oxy— acetylene  mixture.  The  initial  temperature  and  pressure  of  the  mixture 

e 

were  taken  to  be  300  X  and  100  torr  (0.1316  atm)  to  correspond  to  the  initial 
conditions  in  the  experiments  of  Xnystautas  and  Lee6.  As  a  first 


approximation  the  tine  duration  necessary  for  successful  initiation  was 
assumed  to  be  equal  to  the  chemical  induction  time  of  the  mixture 
corresponding  to  the  conditions  at  the  piston  surface.  The  induction  tine 
data  used  were  those  obtained  by  Edwards  et  al. 14  for  an  acetylene-oxygen- 
nitrogen  (2:5:4)  mixture  and  are  given  by: 

i*g  (t[o2]  )  -  -9.41  ±  0.2  +  (6) 

where  t  is  the  induction  tine  in  seconds#  [02]  is  the  concentration  in 
mol /liter,  and  T  is  the  temperature  in  thousands  of  degrees  K.  Three  dif¬ 
ferent  power-energy  relations  obtained  from  the  theoretical  model  are  shown 
in  Figure  1*  Curve  A  was  obtained  by  using  the  smallest  value  of  the  induc¬ 
tion  tins  given  by  Sq.  (6),  that  is#  by  choosing  the  negative  signs.  Curve  B 
was  obtained  by  using  the  mean  values  and  curve  C  by  using  the  largest  value 
of  the  induction  tine  (by  choosing  the  positive  signs).  The  arrows  on  curve 
C  indicate  the  direction  of  increasing  Mach  number.  First#  we  note  that  each 
curve  has  a  minimum  power  and  a  mini  sum  energy.  We  also  observe  that  as  the 
Mach  number  decreases  below  the  Mach  ncmber  corresponding  to  the  min llsnm 
power#  both  the  average  source  power  and  the  source  energy  increase.  How¬ 
ever,  when  the  Mach  number  increases  above  the  Mach  n ember  corresponding  to 
the  minimum  power,  the  energy  first  decreases  to  the  minimum  energy  and  then 
increases  again.  All  three  curves  exhibit  these  sane  qualitative  trends. 

The  shape  of  these  curves  can  be  explained  in  the  following  manner.  As 
the  Need:  number  of  the  shock  wave  decreases#  the  pressure  and  the  temperature 
behind  it  decrease.  This  decrease  also  results  in  a  decrease  of  the  pressure 
and  velocity  at  the  piston  surface.  This  would  tend  to  decrease  both  the 


power  and  the  energy  since,  as  seen  in  Sgs.  (1,2), 

p  **  p  u^t  (7)  and 

P  P 

E  -  p  u2^  .  (8) 

P  P 

nils  tendency  is,  however,  opposed  by  the  tendency  of  the  induction  time  to 
increase  with  decreases  in  the  pressure  and  the  tenperature.  For  low  Mach 
nmbers,  (i.e*,  lew  temperatures  behind  the  shock)  a  small  decrease  in  the 
Mach  number  of  the  shock  wave  leads  to  a  large  increase  in  the  induction 
time*  The  shape  of  the  curves  in  Figure  1  implies  that  this  increase  in 
induction  time  is  more  than  sufficient  to  compensate  for  the  decrease  in  the 
pressure  and  the  velocity  for  Mach  nvmbers  below  that  corresponding  to  the 
minimum  power,  Therefore  both  the  power  and  the  energy  increase  with 
decreasing  Mach  nmriber.  Since  the  energy  is  proportional  to  the  product  of 
the  power  and  the  induction  time,  (Sgs.  (7,8))  the  energy  increases  faster 
with  induction  time  than  the  power  does*  As  the  Mach  number  increases  above 
that  corresponding  to  the  minimi*  power,  the  increase  in  the  pressure  and 
velocity  is  larger  than  the  decrease  in  the  Induction  time.  Therefore  the 
power  increases.  Bowever,  for  a  certin  range  of  Mach  nosbers,  the  increase 
in  the  pressure  and  velocity  is  not  sufficient  to  compensate  for  the  decrease 
in  the  square  of  the  induction  time.  Therefore  the  energy  decreases  until  it 
attains  a  minimum  value,  even  though  the  power  increases.  Finally,  for  Mach 
nxsnbers  above  that  corresponding  to  the  minimum  energy,  the  increase  in  the 
pressure  and  velocity  are  easily  able  to  overcome  the  decrease  in  the  indue** 
tion  time  with  increasing  Mach  nusfeer  and  both  the  power  and  the  energy 
increase.  This  occurs  because  the  rate  of  decrease  of  the  induction  time 


with  temperature  is  small  for  high  temperatures  (i.e.,  high  Mach  numbers) 
according  to  Eq.  (6). 

The  power-energy  curve  obtained  using  data  from  the  spark  ignition 
eiqseriments  ®  of  Knystautas  and  Lee  has  also  been  included  in  Figure  1  as 
curve  D.  This  curve  exhibits  the  same  qualitative  trends  as  those  of  the 
theoretical  curves  discussed  above.  However,  we  observe  that  the  values  of 
the  minimum  power  and  the  ml  nimum  energy  from  the  four  curves  are  very 
different  from  each  other.  Die  differences  in  the  values  of  these  parameters 
from  the  three  "theoretical"  curves  (A,B,  and  C)  indicate  that  the 
experimental  uncertainties  in  the  values  of  the  induction  times  used  have  a 
significant  effect  on  the  value  of  the  minimum  power  and  the  minimum  energy, 
file  ^  power  varies  from  about  0.3  MW/cm  to  about  1  MW/cm  and  the 

■<«<«—  energy  varies  from  about  0.012  J/cm  to  about  0.1  J/csu  The 
experimentally  determined  minimum  power  (from  curve  D)  is  about  0.13  MW/cm, 
which  is  lower  than  the  calculated  values,  and  the  minimum  energy  is  about 
0.1  J/cm,  which  is  at  the  top  of  the  range  of  calculated  values. 

file  quantitative  differences  betwsen  the  exprimental  and  theoretical 
values  could  be  due  to  a  variety  of  factors.  For  example,  there  may  be  dif¬ 
ferences  in  the  shape  of  the  detonation  wave  in  the  experiments  and  the  cal¬ 
culations.  filers  are  also  uncertainties  in  the  induction  time  data.  We  have 
used  the  induction  time  given  by  Bq.  (6)  for  a  range  of  t superstores  far 
greater  than  that  over  which  it  was  determined.  As  seen  from  Figure  1  and 
the  discussion  above,  uncertainties  in  the  induction  time  have  a  significant 
effeet  on  the  power-energy  relations.  Furthermore  some  of  the  approximation* 
mad*  in  the  theoretical  model,  such  as  the  assonption  that  the  similarity 
approximation  is  valid  or  that  the  appropriate  time  to  be  used  in  lq.  (1,2) 


is  the  induction  tine,  are  also  debatable.  These  factors  and  other  limit* 
ations  of  the  model  have  been  studied  and  discussed  in  more  detail  by 
Kailasanath  and  Oran 11 . 

The  derived  power-energy  relation  for  the  initiation  of  planar  detona¬ 
tions  in  the  same  oxy -acetylene  mixture  is  shown  in  Figure  2.  In  this 
figure,  we  also  show  the  shock  tube  data  of  Dabora7  on  the  direct  initiation 
of  detonations  in  a  stoichiometric  hydrogen-air  mixture.  Both  curves  exhibit 
the  same  qualitative  behavior.  Alike  the  cylindrical  case,  each  value  of 
the  power  corresponds  to  an  unique  value  of  energy.  The  direction  of 
increasing  shock  strength  (as  determined  by  the  Mach  number)  is  also  shown  in 
Figure  2.  In  the  planar  case,  we  see  that  as  the  Mach  number  decreases  the 
power  always  decreases.  As  noted  earlier,  in  the  cylindrical  ease,  as  the 
Mach  number  decreases,  the  power  decreases  only  up  to  the  minimum  power. 

Then  the  power  increases  with  a  decrease  in  the  Mach  number  of  the  shock 
wave.  Therefore,  the  qualitative  difference  in  the  ejperimental  data  of 
Xhystautas  and  lee  (shown  in  Figure  1)  and  Dabora  (shown  in  Figure  2)  are  due 
to  the  difference  in  the  geometry  of  the  two  detonations. 

In  Figure  2,  we  also  observe  that  as  the  Mach  nonber  decreases,  we  need 
more  and  sore  energy  to  Initiate  a  detonation.  The  trend  of  the  curves 
indicates  that  there  is  a  minimum  Mach  number  below  which  a  detonation  will 
not  occur  (i.e. ,  would  require  an  infinite  amount  of  energy).  The  value  of 
the  power  corresponding  to  this  mlnlncm  Mach  number  is  the  minimum  power. 

This  agrees  with  the  observation  ande  by  lee  et  al.  6  that  the  requirement  for 
a  minimum  value  of  the  source  power  indicates  that  the  source  must  be  capable 
of  generating  a  shook  wave  of  a  certain  minimum  Mach  ntaber.  However,  we 
observe  from  Figure  1  that  for  the  case  of  cylindrical  detonations,  the 
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■ini—  power  does  not  correspond  to  the  shock  wave  of  min  liman  Mach  number 


In  the  cylindrical  case/  it  is  possible  to  initiate  a  detonation  with  a  shock 
wave  of  lower  Mach  umber  than  that  corresponding  to  the  minimum  power.  Such 
a  shock  will  have  to  be  maintained  for  a  longer  time  than  the  shock  corres¬ 
ponding  to  the  minimum  power  and  hence  will  require  a  larger  amount  of 
energy. 

The  power-energy  curve  for  the  initiation  of  spherical  detonations  is 
similar  to  the  curve  for  the  cylindrical  case.  However,  for  the  ease  of 
spherical  detonations,  the  power, 


P  -  p  u3  t2 
P  P 

but  we  still  have  the  energy. 


(9) 


E  *  P  t  (10) 

Since  the  power  and  energy  are  proportional  to  higher  powers  of  the  time,  t, 
uncertainties  in  t  will  have  a  greater  effect  on  the  value  of  the  minimum 
power  and  the  minimum  energy.  Further  work  is  being  carried  out  currently  to 
study  the  initiation  of  spherical  detonations  in  hydrogen-air  mixtures  and  to 
oospare  it  with  ei^erimental  data. 


Conclusions 

the  results  discussed  above  show  that  though  the  single  theoretical 
model  has  significant  limitations,  it  can  still  be  used  to  explain  the  quali¬ 
tative  differences  in  the  power-energy  relations  obtained  from  different 
experimental  arrangements.  The  results  from  the  model  also  shew  that  the 
minimum  power  and  the  minimum  detonation  energy  are  very  sensitive  to  uncer¬ 
tainties  in  the  ignition  delay  time  used,  toother  interesting  conclusion 


from  the  model  is  that  the  minimum  power  requirement  corresponds  to  a  shock 
of  minimum  Mach  number  only  in  the  case  of  planar  detonations.  Currently, 
further  applications  of  the  model  and  extensions  to  the  study  of  condensed 
phase  detonations  are  being  considered. 
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Figure  1.  Power-energy  relations  for  the  initiation  of  cylindrical 
detonations  in  an  ecetylen e-oxygen -nitrogen  mixture 
(2:5:4)  at  1  ata  and  300  K.  Curve  D  is  obtained  from 
spark  ignition  cxperlaents  (6).  Curves  A,  B,  and  C  are 
explained  in  the  text.  The  arrows  on  Curve  C  indicate  the 
direction  of  increasing  Mach  number. 
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Power-energy  relations  for  the  initiation  of  planar 
detonations.  The  x'a  are  data  obtained  froa  shock 
tube  experiaenta  (7). 


Numerical  Simulations  of  Weak  Ignition  Behind  Incident  Shock  Waves 
and  the  Transition  to  Detonation 

K.  Kailasanath*  and  E.S.  Oran 
Laboratory  for  Computational  Fhysics 
Naval  Research  Laboratory 
Washington,  DC  20375 


Introduction 

Shock  induced  ignition  in  gaseous  hydrogen-air  mixtures  may  occur  in  one 
of  two  distinct  modes,  depending  on  the  thermodynamic  state  in  the  shocked 
material.  At  lover  tesperatures  the  ignition  is  veak,  or  mild,  with  the 
gradual  development  of  the  gas  dynamic  explosion;  at  higher  tenperatures  it 
is  strong,  or  sharp,  with  a  practically  instantaneous  appearance  of  a 
secondary  shock  induced  by  the  explosive  reaction.  Furthermore  at  lov 
tenperatures  the  formation  of  distinct  flame  kernels  appears  to  be  an 
essential  precursor  of  ignition.  Meyer  and  Oppenheim  [l )  point  out  that 
the  intrinsically  turbulent  flow  field  behind  a  reflected  shock  in  a  shock 
tube  results  in  a  nonuniform  temperature  distribution  which  creates  distinct 
reaction  centers  and  leads  to  weak  ignition.  Although  the  nonuniformity  of 
tesperature  caused  by  turbulence  is  one  mechanism  which  gives  rise  to 
reaction  centers,  the  nonsteadiness  in  the  velocity  of  the  causal  shock  can 
also  produce  them.  This  latter  effect  has  been  convincingly  demonstrated  by 
9trehlov  et  al.  [2]  in  their  studies  of  shock  propagation  in  a  slowly 
converging  channel.  It  haB  also  been  shown  that  when  an  incident  shock  in  a 
uniform  shock  tube  begins  to  accelerate,  hot  spots  occur  before  the 
transition  to  detonation  [3,M* 

We  present  results  from  time  dependent  numerical  sinulations  used  to 
elucidate  some  of  the  details  of  veak  ignition  behind  incident  shock  waves 
and  the  subsequent  transition  to  detonation.  These  simulations  were 
performed  using  a  one-dimensional  model  {5,6]  which  combines  a  description 
of  the  fluid  dynamics  and  detailed  chemical  kinetics.  With  such  detailed 
one-dimensional  numerical  sisulations  one  hopes  to  isolate  and  study  some  of 
the  various  phenomena  occuring  in  actual  nulti-dimensional  systems. 

The  Numerical  Model 


The  one-klimensional  reactive  shock  model  ( 5 ,6]  used  to  perform  the 
calculations  described  below  solves  the  time-dependent  conservation 
equations  |7]  for  mass,  momentum  and  ener^  coupled  to  the  equations 
describing  the  chemical  kinetics.  The  model  uses  an  explicit,  Eulerian 
finite  difference  formilation  with  a  sliding  rezone  capability  to  provide 
resolution  around  moving  gradients.  The  solutions  of  the  equations 
describing  the  fluid  dynamics  and  the  chemistry  of  the  problem  are  coupled 
using  time-step  splitting  techniques  |6]. 

The  convective  transport  terms  in  the  conservation  equations  are  solved 
— using  one  variant  of  the  Flux-Corrected  Transport  (FCT)  method'  (8] .  This  IF 
~  a  conservative,  monotonic  algorithm  with  fourth-order  phase  accuracy  and 
does  not  require  artificial  viscosity  to  stabilize  shocks.  The  ordinary 
differential  equations  describing  the  chemical  kinetics  are  solved  using 

) 
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VSALM,  a  vectorized  version  of  the  selected  asymptotic  integration  method 
enplcyed  in  CHEMEQ  [9].  The  chemical  kinetics  rate  scheme  used  consists  of 
about  fifty  rates  relating  the  species  E2,  02,  H,  0,  OH,  H02,  H20  and  H202 
and  has  been  extensively  tested  against  experimental  data  [10 ,11] . 

Results  and  Discussion 


In  the  shock  tube  s initiations,  the  pressure  across  the  diaphragm  and  the 
driven  gas  mixture  are  chosen  so  that  the  thermodynamic  state  behind  the 
incident  shock  is  in  the  weak  ignition  regime  [12].  Figure  1  is  a  position¬ 
time  diagram  of  the  events  occurring  in  the  shock  tube  simulation.  The 
trajectory  of  the  shock  front  is  labelled  S  and  that  of  the  contact  surface 
is  labelled  CS.  Five  different  regimes  have  also  been  identified  on  the 
diagram.  They  are  (a)  pre-ignition  regime,  (b)  quasi-steady  shock-reaction 
couplex,  (c)  formation  of  reactive  centers,  (d)  hot  spot  formation  leading 
to  an  overdriven  detonation  and  (e)  detonation  relaxation.  Similar  regions 
have  been  Identified  by  Edvards  et  al.  [1]  in  their  shock  tube  experiments. 
These  regions  are  examined  in  detail  below. 

Initially,  the  shock  travels  at  a  nearly  constant  velocity  (1.1»  x  10s 
em/s)  into  the  hydrogen-air  mixture,  raising  its  temperature  and  pressure  to 
a  near-constant  value.  Gradually  reactions  begin  to  occur  in  the  shock 
heated  gam.  Small  pressure  disturbances  due  to  energy  release  travel 
forward  at  a  velocity  of  1.8  X  10s  cm/s  (in  the  laboratory  frame  of 
reference),  which  is  the  sum  of  the  sonic  and  particle  velocities  behind  the 
shock.  When  they  reach  the  shock  they  accelerate  the  shock  slightly.  Cbee 
the  shock  is  accelerated,  the  tenperature  of  the  gases  now  passing  through 
the  shock  is  higher  than  the  raised  teaperature  created  by  the  original 
shock.  In  the  veak-ignltion  regime,  such  an  increase  in  temperature  can 
result  in  a  significant  reduction  in  the  induction  time  [12,  13].  The 
numerical  simulations  shew  how  this  leads  to  the  formation  of  reactive 
centers  (RC)  in  the  nevly  shocked  material.  At  these  reactive  centers, 
reaction  progresses  at  a  ouch  move  rapid  rate  than  in  the  previously  shocked 
material.  The  development  of  a  hot  spot  due  to  energy  release  at  one  of  the 
reactive  centers  can  be  seen  in  the  simulations.  The  presence  of  such  hot 
spots  has  been  observed  earlier  in  shock  tube  experiments  [2,M» 

Energy  release  at  a  hot  spot  causes  a  pair  of  flamelets  or  reaction 
vaves  (RW),  one  propagating  ahead  into  the  shock  heated  gas  mixture  and  the 
other  propagating  back  towards  the  contact  surface.  These  flamelets 
initially  propagate  at  a  subsonic  speed  with  respect  to  the  fluid.  Such 
flamelets  have  been  observed  by  Bazhenova  and  Soloukhin  (3l  in  their 
Incident  shock  tube  experiments.  Iher&  release  behind  these  flamelets 
causes  pressure  vaves  vhich  accelerate  the  flamelets  into  detonation  type 
vaves.  The  reaction  vave  moving  backward  against  the  fluid  accelerates 
slowly  and  travels  at  nearly  the  sonic  speed  till  it  reaches  the  contact 
surface.  However  the  forward  moving  reaction  vave  transitions  into  a 
strong  detonation  vave  even  before  it  reaches  the  incident  shock  vave.  This 
agrees  with  the  observation  made  by  Bazhenova  and  Soloukhin  [3]  that  the 
merging  of  the  flame  front  with  the  incident  shock  vave  is  not  a  necessary 
condition  for  the  detonation  vave  formation. 
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The  numerical  simulations  presented  in  this  paper  shew  that  the  one- 
dimensional  reactive  flow  model  with  detailed  chemical  kinetics  can  be  used 
to  elucidate  some  of  the  details  of  weak  ignition  behind  incident  shocks  and 
the  subsequent  transition  to  detonation.  Die  predictions  of  the  model  are 
in  qualitative  agreement  with  experimental  observations.  However,  the  model 
does  not  include  suit i -dimensional  phenomena  such  as  turbulence  and  boundary 
layer  growth  which  play  an  important  part  in  any  quantitative  study  of  the 
transition  to  detonation.  Currently  a  two-dimensional  reactive  shock  model 
exists  but  it  uses  a  parameterized  model  for  ener^1  release  ll4,15).  The 
model  is  now  being  extended  to  include  a  detailed  chemical  kinetic  scheme. 
Calculations  such  as  shown  here  performed  with  this  model  would  show  the 
effects  of  transverse  waves  and  boundary  layers  on  the  transition  to 
detonation. 

This  work  has  been  supported  hy  the  Office  of  Naval  Research  through  the 
Naval  Research  Laboratory. 
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Figure  1.  A  position-tin*  diagram  of  the  main  events  occurring 
in  the  shock  tube  simulation.  The  symbols  are 
explained  in  the  text. 
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ABSTRACT 


Time-dependent  numerical  simulations  are  used  to  elucidate  some  of  the 


details  of  weak  ignition  behind  incident  shocks  and  the  subsequent  transition 
to  detonation.  It  is  shown  that  a  small  amount  of  energy  released  in  the 
shocked  region  can  be  the  origin  of  pressure  waves  which  accelerate  the  shock 
front*  The  simulations  presented  here  show  how  this  leads  to  the  formation 
of  reactive  centers.  The  for  nation  of  a  hot  spot  due  to  energy  release  at 
one  of  the  reactive  centers  and  the  subsequent  development  of  a  pair  of 
flamelets  from  the  hot  spot  are  studied  using  the  numerical  simulations.  The 
results  of  the  simulations  are  also  compared  to  experimental  observations. 
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I.  IBTRODUCTIOR 


Shock-induced  ignition  in  gaseous  hydrogen-air  mixtures  may  occur  in  one 
of  tvo  distinct  modes,  depending  on  the  theraotynamic  state  in  the  shocked 
material.  At  lower  temperatures  the  ignition  is  weak,  or  mild,  with  the 
gradual  development  of  the  gas  tynamic  explosion.  At  higher  temperatures  it 
is  strong,  or  sharp,  with  an  abrupt  appearance  of  a  secondary  shock  induced 
by  the  exploslre  reaction.  Furthermore  at  low  temperatures  the  formation  of 
distinct  flame  kernels  appears  to  be  an  essential  precursor  of  ignition. 

Meyer  and  Cppenheim  (1971a)  point  out  that  the  intrinsically  turbulent  flow 
field  behind  a  reflected  shock  in  a  shock  tube  results  in  a  nonuniform 
tesperature  distribution  which  creates  distinct  reaction  centers  and  leads  to 
weak  ignition.  Although  the  nonuni f oradty  of  tesperature  caused  by 
turbulence  is  one  mechanism  which  gives  rise  to  reaction  centers,  the 
nonsteadiness  in  the  velocity  of  the  causal  shock  can  also  produce  them. 

This  latter  effect  has  been  convincingly  demonstrated  by  Strehlcw  et  al. 
(1967)  in  their  studies  of  shock  propagation  in  a  slowly  converging  channel. 
It  has  also  been  shown  that  when  an  incident  shock  in  a  uniform  shock  tube 
begins  to  accelerate,  hot  spots  occur  before  the  transition  to  detonation 
(Bazhenova  and  Soloukhin,  1959;  Edwards  et  al. ,  198l). 

In  this  paper  the  reactive  flow  behind  an  incident  shock  wave  in  a 
hydrogen-air  mixture  is  simulated  using  a  one-dimensional,  time -dependent 
numerical  model  which  combines  a  description  of  the  fluid  dynamics  and 
detailed  chemical  kinetics.  In  the  slnulations,  the  pressure  ratio  across 
the  diaphragm  and  the  driven  gas  mixture  are  chosen  so  that  the  thermodynamic 
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state  behind  the  incident  shock  is  in  the  weak  ignition  regime  (Oran  and 
Boris,  1981b).  It  is  shown  that  a  small  amount  of  energy  released  in  the 
shocked  gas  (which  might  occur  due  to  density,  temperature  or  stoichiometric 
fluctuations)  can  be  the  origin  of  pressure  waves  which  accelerate  the  shock 
front.  Once  the  shock  is  accelerated,  the  temperature  of  the  gases  now 
passing  through  the  shock  is  higher  than  the  raised  temperature  created  by 
the  original  shock.  In  the  weak  ignition  regime,  such  an  increased 
teaperature  can  result  in  a  significant  reduction  in  the  induction  time 
(Meyer  and  Oppenheim.  1971b;  Oran  and  Boris,  1981b).  The  simulations 
presented  here  show  how  this  leads  to  the  formation  of  reactive  centers  in 
the  newly  shocked  material  where  reaction  progresses  at  a  more  rapid  rate 
than  in  the  previously  shocked  material.  The  formation  of  a  hot  spot  due  to 
energy  release  at  one  of  the  reactive  centers  and  the  subsequent  development 
of  a  pair  of  flamelets  or  reaction  waves  from  the  hot  spot  are  studied  using 
the  numerical  simulations.  The  results  of  the  simulations  have  also  been 
compared  to  experimental  observations  (Bazhenova  and  Soloukhin,  1959;  Urtiew 
and  Oppenheim,  1966,  1967;  Strehlow  et  al. ,  1967;  Edwards  et  al. ,  19SL/. 

II.  THE  NUMERICAL  MODEL 

The  one-dimensional  reactive  shock  model  (Oran  et  al. ,  1979;  Oran  and 
Boris,  1981a)  used  to  perform  the  calculations  described  below  solves  the 
time-dependent  conservation  equations  (Williams,  1965;  Oran  and  Boris,  1981a) 
for  mass,  momentum  and  energy  coupled  to  the  equations  describing  the 
chemical  kinetics.  The  model  uses  an  explicit,  Eulerlan  finite  difference 
formulation  with  a  sliding  rezone  capability  to  provide  resolution  around 
moving  gradients.  The  solutions  of  the  equations  describing  the  fluid 
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dynamics  and  the  chemistry  of  the  problem  are  coupled  using  time-step 
splitting  techniques  (Oran  and  Boris,  1981a). 

The  convective  transport  terms  in  the  conservation  equations  are  solved 
using  one  variant  of  the  Flux-Corrected  Transport  (FCT)  method  (Boris  and 
Book,  1976;  Boris,  1976).  This  is  a  conservative,  monotonic  algorithm  vith 
fourth-order  phase  accuracy  and  does  cot  require  artificial  viscosity  to 
stabilize  shocks.  The  Flux  Correction  procedure  itself  ensures  that  the 
shocks  are  one  or  two  zones  vide  and  have  maximal  resolution.  The  ordinary 
differential  equations  describing  the  chemical  kinetics  are  solved  using 
VS AIM,  a  vectorized  version  of  the  selected  asymptotic  integration  method 
employed  in  CH1MEQ  (Young  and  Boris,  1977;  Young,  I9B0).  This  algorithm 
identifies  the  stiff  equations  for  treatment  vith  a  stiffly  stable  method. 

The  remaining  equations  are  solved  vith  a  standard  classical  method.  The 
algorithm  has  been  specially  optimized  for  use  in  conjunction  vith  fluid 
dynamic  models. 

The  chemical  kinetics  rate  scheme  used  is  given  in  Table  I.  It  consists 
of  about  fifty  rates  relating  the  species  H2,  02*  H,  0,  OH,  H02,  H20  and 
and  has  been  extensively  tested  against  experimental  data  (Oran  et  al. ,  19&L, 
Burks  and  Oran,  1980).  Burks  and  Oran  (1980)  showed  that  the  results 
computed  vith  the  scheme  compared  very  veil  vith  experimentally  observed 
Induction  times,  second  explosion  limits  and  the  temporal  behavior  of 
reactive  species.  Oran  et  al.  (19SO  have  shown  that  the  scheme  gives  good 
results  when  coupled  vith  a  fluid  dynamic  model  in  the  simulation  of  the 
conditions  behind  a  reflected  shock.  Beats  of  formation  and  enthalpies  have 
been  taken  from  the  JAHAF  tables  (1971). 


C-6 


For  the- calculations  performed  in  this  paper,  the  timescales  under 
consideration  are  short  and  therefore  the  diffusive  transport  processes, 
thermal  conduction  and  molecular  diffusion,  have  negligible  effect.  The 
effects  of  these  processes  have  not  been  considered  although  they  are  part  of 
the  general  numerical  model.  Although  the  geometry  may  be  either  cartesian, 
cylindrical,  spherical  or  some  generalized  co-ordinate,  the  simulations 
presented  below  are  in  cartesian  geometry. 

The  detailed  simulations  discussed  in  this  paper  require  that  ve  model 
relatively  long  systems  (on  the  order  of  meters)  while  we  simultaneously 
maintain  high  accuracy  around  steep  gradients  such  as  the  shock  front  and  the 
contact  surface.  The  rather  sophisticated  adaptive  gridding  method  developed 
for  this  purpose  is  shown  schematically  in  Figure  1.  There  are  two  finely 
gridded  regions:  One  surrounding  the  shock  wave  and  the  other  surrounding 
the  contact  surface.  The  fine-zoned  region  around  the  contact  surface  moves 
with  the  contact  surface  at  the  fluid  velocity,  and  so  this  part  of  the 
calculations  is  essentially  Iagrangian.  The  region  around  the  shock  front 
moves  with  the  front.  Jhch  of  these  finely  gridded  regions  may  have  a 
different  minimum  computational  cell  size.  The  computational  cells  in  the 
regions  ahead  of  the  shock  wave  and  behind  the  contact  surface  change 
exponentially  in  size  from  the  smallest  near  the  shock  wave  or  the  contact 
surface  to  the  largest  at  the  walls.  Care  is  taken  that  the  transition  in 
the  cell  sizes  is  smooth.  Fbr  the  results  presented  in  this  paper  a  total  of 
200  computational  cells  are  used  to  describe  the  shock  tube  and  the  cell 
sizes  varied  from  0.1  cm  around  the  shock  to  over  50  cm  near  the  shock  tube 
end-walls . 
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III.  WEAK  IGNITION  BEHIND  AN  INCIDENT  SHOCK 


The  numerical  model  described  in  the  preceding  section  was  used  to  study- 
weak  ignition  behind  planar  incident  shock  waves.  The  system  parameters,  the 
initial  conditions  and  the  temperature  and  pressure  behind  the  incident  shock 
have  all  been  summarized  in  Table  II.  The  chemical  induction  time,  for  the 
conditions  described  by  the  incident  shock,  is  about  2000  -s.  However,  since 
the  thermodynamic  state  behind  the  shock  is  in  the  weak  ignition  regime, 
ignition  may  occur  much  earlier  than  2000  "s  due  to  temperature,  pressure,  or 
density  fluctuations. 

Figure  2  is  a  position-time  diagram  of  the  events  occurring  in  the  shock 
tube  simulation.  The  trajectory  of  the  shock  front  is  labelled  S  and  that  of 
the  contact  surface  is  labelled  CS.  Except  for  small  variations  (which  are 
examined  in  detail  below)  the  shock  travels  at  a  nearly  constant  velocity, 

1.1  x  10s  cm/s,  until  the  reaction  wave  formed  between  the  contact  surface 
and  the  shock  front  reaches  it.  At  this  time  the  velocity  rises  quickly  to 
3.21  x  ics  cm/s.  It  then  gradually  decreases  towards  the  Chapman-Jouguet 
detonation  velocity.  Five  different  regimes  have  also  been  identified  on  the 
diagram.  They  are  (a)  pre-ignition  regime,  (b)  quasi-steady  shock-reaction 
complex,  (c)  formation  of  reactive  centers,  (d)  hot  spot  formation  leading  to 
an  overdriven  detonation  and  (e)  detonation  relaxation.  Similar  regions  have 
been  identified  by  Edwards  et  al.  (1961)  in  their  shock  tube  experiments. 
These  regions  are  examined  in  detail  below. 

As  the  shock  travels  at  a  nearly  constant  velocity  into  the  hydrogen-air 
mixture,  the  temperature  and  pressure  of  the  mixture  are  raised  to  a  near¬ 
constant  value.  Reactions  in  the  shock  heated  gas  first  occur  near  the 


contact  surface  since  the  temperature  has  been  high  for  the  longest  time 
here,  as  has  been  observed  by  Urtiew  and  Oppenheim  (1967).  In  Figure  3,  the 
spatial  variation  of  the  temperature  and  the  OE  mole  fraction  between  the 
shock  and  the  contact  surface  are  shown  at  135  ys  after  the  bursting  of  the 
diaphragm.  The  OH  mole  fraction  attains  a  maximum  value  near  the  contact 
surface  and  decreases  across  the  system  towards  the  shock  front.  This  is 
because  the  gas  mixture  at  different  locations  has  been  at  the  higher 
temperature  for  different  durations. 

Formation  of  Reactive  Centers  and  Hot  Spots 

Small  pressure  disturbances  occur  due  to  the  energy  released  in  the 
reactions  near  the  contact  surface.  These  pressure  disturbances  travel 
forward  at  a  velocity  (in  the  laboratory  frame  of  reference)  of  1.8  x  105 
cm/s,  which  Is  the  sum  of  the  sonic  and  particle  velocities  behind  the  shock. 
When  they  reach  the  shock  they  accelerate  the  shock  slightly  resulting  in  a 
temperature  increase  behind  the  shock.  This  can  be  observed  in  Figure  U 
where  the  spatial  distribution  of  the  temperature  between  the  shock  and  the 
contact  surface  is  shown  at  four  different  times.  The  nonuniform  temperature 
distributions  can  be  explained  by  the  following  sequence  of  events:  pressure 
disturbances  originating  at  different  times  near  the  contact  surface  reach 
the  incident  shock  at  different  times  and  each  successive  pressure 
disturbance  meets  a  shock  of  slightly  different  strength  because  the  shock 
has  already  been  accelerated  slightly  due  to  previous  pressure  disturbances. 
We  know  that  in  the  weak  ignition  regime,  the  induction  time  is  vezy 
sensitive  to  perturbations  in  teiqperature  and  pressure  (Oran  and  Boris, 
19&lb).  This  sensitivity,  in  fact,  is  what  produces  the  two  peaks  in  OH  mole 


fraction  distribution  (Figure  5A)  at  3^8  us  after  the  bursting  of  the 
diaphragm.  The  first  peak  is  due  to  the  development  of  reactions  near  the 
contact  surface  and  occurs  in  the  same  reactive  center  (see  Figure  2)  vhich 
was  observed  earlier.  Nov  a  second  reactive  center  has  formed  closer  to  the 
shock  front  due  to  the  higher  temperature  created  by  the  acceleration  of  the 
shock.  Again  ve  observe  that  the  reactive  center  occurs  in  that  part  of  the 
shock  heated  gas  vhich  vas  at  the  elevated  temperature' for  the  longest  period 
of  time.  Ahead  of  the  second  reactive  center  the  rapid  decrease  in  the 
radical  mole  fraction  is  halted  by  a  further  increase  in  the  temperature 
behind  the  shock  and  a  third  reactive  center  is  forming.  By  UU5  vs  (Figure 
5B),  a  third  reactive  center  has  developed  and  the  radical  mole  fraction  is 
rising  most  rapidly  at  this  center.  The  temperature  around  this  reactive 
center  is  significantly  (about  Uo°  K)  greater  than  that  behind  the  initial 
incident  shock  vave.  In  the  veak  ignition  regime,  such  changes  in 
temperature  result  in  a  substantial  reduction  in  the  induction  time. 

Therefore  reaction  progresses  at  this  reactive  center  at  a  much  more  rapid 
rate  than  at  any  other  location  in  the  system.  This  results  in  a  "hot  spot" 
in  the  system.  The  development  of  this  hot  spot  can  be  seen  in  Figure  6 
where  the  temperature -distribution  between  the  contact  surface  and  the  shock 
front  has  been  shown  for  four  different  times.  The  temperature  increases  ty 
more  than  200°  in  30  us. 

Transition  to  Detonation 

The  events  occurring  after  the  hot  spot  formation  can  be  seen  in  Figure 
7  where  a  segment  of  the  position-time  diagram  (Figure  2)  is  presented  in 
greater  detail.  The  hot  spot  (marked  "A"  in  the  figure)  travels  with  the 


fluid  and  continuously  releases  ener©-  into  it.  This  results  in  "flamelets" 
or  "reaction  waves"  ("B")  which  travel  with  respect  to  the  fluid.  One 
reaction  wave  moves  forward  with  the  fluid  while  another  moves  against  the 
fluid  towards  the  driver  section.  These  reaction  waves  initially  move 
subsonically  with  respect  to  the  fluid  hut  are  soon  accelerated  into  steep 
detonation  type  waves  which  move  supersonically  with  respect  to  the  fluid. 
When  the  forward  moving  reaction  wave  reaches  the  shock  front,  the  shock 
velocity  abruptly  increases  to  a  high  value.  The  shock-reaction  conplex  then 
moves  as  a  strong  overdriven  detonation  wave  which  decelerates  towards  the 
Chapman-Jouguet  detonation  velocity. 

The  temperature,  pressure  and  velocity  distributions  across  the  system 
at  a  particular  time  after  the  hot  spot  has  formed  are  shown  in  Figure  8a. 

By  this  time  the  energy  release  has  resulted  in  a  noticeable  pressure  rise. 
This  occurs  because  the  energy  release  is  occurring  at  nearly  constant  volume 
conditions.  The  energy  release  at  the  hot  spot  causes  a  series  of  minute 
pressure  pulses  to  propagate  both  forward  and  backward,  each  a  little 
stronger  than  the  previous  one.  A  series  of  pressure  pulses  are  produced 
since  the  energy  release  occurs  over  a  period  of  time  and  is  determined  at 
each  time  from  the  detailed  chemical  kinetic  interactions  among  the  various 
species.  These  pressure  pulses  coalesce  to  form  steepening  pressure  waves 
propagating  into  the  shocked  mixture  as  seen  in  Figure  8b.  The  time  history 
of  these  reaction  waves  has  been  shown  in  a  series  of  figures  (Fig.  8b  - 
Fig.  10b)  which  cover  the  time  period  between  609  us  and  688  us. 

Let  ub  first  look  at  the  development  of  the  reaction  wave  which  moves 
forward  into  the  shock-heated  gas  mixture.  The  velocity  of  the  forward 
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moving  vave,  the  fluid  velocities  on  both  sides  of  the  wave  and  the  speed  of 
sound  in  the  gas  mixtures  on  both  sides  of  the  wave  are  given  in  Thble  III 
for  a  series  of  times.  At  609  us,  the  forward  moving  reaction  wave  is 
supersonic.  Die  velocity  of  the  fluid  with  respect  to  the  wave  decreases 
across  the  vave  but  is  supersonic  on  either  side.  Furthermore  the  pressure 
rise  across  the  vave  is  just  over  a  factor  of  two,'  which  is  moderate.  Dins 
the  reaction  vave  behaves  like  a  weak  detonation  vave  at  this  time.  later, 
as  seen  in  Figure  9s,  the  pressure  rise  across  the  vave  has  increased  and  the 
wave  is  also  travelling  faster.  Die  velocity  of  the  fluid  with  respect  to 
the  wave  is  still  supersonic  on  both  sides  of  the  wave.  The  pressure  rise 
across  the  vave  continues  to  increase  and  at  the  time  corresponding  to  Figure 
9b,  the  fluid  velocity  behind  the  wave  is  nearly  sonic.  Die  weak  detonation 
seems  to  be  transitioning  into  a  strong  detonation.  It  does  so  later  (Figure 
9c)  when  the  fluid  velocity  changes  from  supersonic  to  subsonic  across  the 
reaction  vave.  Die  pressure  rise  across  the  vave  is  also  larger  now.  Die 
observed  acceleration  of  the  forward  moving  vave  into  a  strong  detonation  is 
due  to  the  nonlinear  interaction  between  chemical  kinetics  and  fluid 
dynamics.  When  the  forward  moving  vave  moves  into  the  previously  shocked 
material  there  is  a  large  pressure  and  temperature  rise  across  it  since  it  is 
a  strong  compressive  vave.  Diis  increase  in  the  pressure  and  the  temperature 
reduces  the  induction  time  of  the  material  which  crossed  the  wave.  Ihergy 
release  in  this  newly  re-compressed  material  accelerates  the  forward  moving 
vave  further,  and  this  cycle  is  repeated  until  the  forward  moving  wave 
reaches  the  incident  shock  vave  (Figure  10a).  By  this  time  the  pressure 
spike  behind  the  reaction  vave  has  risen  to  6.6  *  106  dynes/cm2.  Because  the 


reaction  wave  is  moving  towards  the  incident  shock  at  a  speed  greater  than 
the  local  speed  of  sound  in  the  shock-heated  gas,  there  is  no  advance  warning 
of  the  over-pressure  region  until  the  main  spike  physically  arrives  at  the 
shock  front.  When  the  reaction  wave  does  coalesce  with  the  shock,  there  is  a 
very  rapid  increase  in  the  shock  speed  (Figure  7)  and  the  shock -react ion  wave 
complex  moves  as  an  overdriven  strong  detonation  wave.  Itee  to  the  high  over¬ 
pressure  that  was  associated  with  the  reaction  front,  the  detonation  wave 
overshoots  the  Ch&pnan-Jouguet  value  by  a  substantial  amount.  However,  the 
overshoot  cannot  be  sustained  by  the  reaction  and  the  overdriven  detonation 
gradually  relaxes  t cwards  a  Chapman- Jouguet  wave.  By  688  us  (Figure  10b)  the 
pressure  spike  behind  the  detonation  vave  has  decreased  to  1*.8  *  10* 
dynes/cm2  and  the  wave  is  moving  faster  than  the  C-J  velocity  by  only  36%. 

In  Figure  10b  we  also  observe  a  small  amplitude  pressure  wave  moving  into  the 
detonation  products.  This  pressure  wave  was  formed  when  the  reaction  wave 
interacted  vith  the  shock  front. 

Let  us  now  look  at  the  reaction  vave  which  moves  backward  towards  the 
driver  section.  From  Figure  7  and  Table  IV  we  see  that  this  wave  moves  at 
nearly  the  sonic  speed.  The  fluid  velocity  vith  respect  to  the  wave  is 
supersonic  on  either  side  of  the  vave.  The  pressure  rise  and  the 
acceleration  of  the  wave  are  slower  since  the  wave  is  moving  against  the 
fluid.  The  pressure  increase  behind  the  wave  is  broader  (Figure  9c)  and 
smaller  than  that  of  the  forward  moving  wave.  It  continues  to  propagate  like 
a  weak  detonation  wave  until  it  interacts  with  the  contact  surface  (Figure 
IQs).  This  interaction  produces  a  pressure  pulse  which  travels  into  the 
helium  driver  gas.  In  Figure  10b  we  see  that  this  pressure  pulse  has 
produced  a  slight  temperature  increase  in  the  helium. 
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IV.  DISCUSSION  AND  CONCLUSIONS 


In  this  paper  the  reactive  flow  behind  an  incident  shock  wave  has  been 
studied  using  a  time -dependent  numerical  model  which  includes  both  detailed 
chemical  kinetics  and  one-dimensional  fluid  dynamics.  The  numerical 
simulations  show  that  the  incident  shock  initially  travels  at  a  steady  speed 
leaving  behind  it  material  in  what  has  been  called  a  "quasi-steady  reaction 
complex"*  From  their  extensive  experimental  studies  on  detonation,  Lee  et 
al.  (1976*  1977)  have  concluded  that  such  regimes  are  universal  predetonation 
phenomena.  This  regime  has  also  been  observed  in  the  incident  shock  tube 
experiments  of  Edvards  et  al.  (1981). 

The  calculations  have  also  shown  the  acceleration  of  the  shock  due  to 
pressure  waves  created  by  energy  release  in  the  shock  heated  gas  mixture,  a 
phenomenon  also  observed  by  Edvards  et  al.  (19&L).  Shock  acceleration  raises 
the  temperature  of  the  gases  passing  through  the  shock  and  in  the  weak 
ignition  regime,  this  results  in  a  significant  reduction  of  the  induction 
time  (Meyer  and  Oppenheim,  1971b;  Oran  and  Boris,  1981b).  This  leads  to  the 
formation  of  reactive  centers  where  reaction  progresses  at  a  more  rapid  rate 
than  in  the  previously  shocked  material.  The  development  of  a  hot  spot  due 
to  energy  release  at  one  of  the  reactive  centers  has  been  shown  in  Figure  6. 
The  presence  of  such  hot  spots  has  been  observed  earlier  in  shock  tube 
experiments  (Strehlov  et  al. ,  1967,  Edvards  et  al. ,  1981). 

Xhergy  release  at  a  hot  spot  causes  a  pair  of  flamelets  or  reaction 
waves,  one  propagating  ahead  into  the  shock  heated  gas  mixture  and  the  other 
propagating  back  towards  the  contact  surface.  These  flamelets  initially 
propagate  at  a  subsonic  speed  with  respect  to  the  fluid.  Such  flamelets  have 
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been  observed  by  Bazhenova  and  Soloukhin  (1959)  in  their  incident  shock  tube 
experiments.  Eher©’’  release  behind  these  flanelets  causes  pressure  waves 
which  accelerate  the  flanelets  into  detonation  type  waves.  The  reaction  wave 
moving  backward  against  the  fluid  accelerates  slowly  and  travels  at  nearly 
the  sonic  speed  till  it  reaches  the  contact  surface.  However  the  forward 
moving  reaction  wave  transitions  into  a  strong  detonation  wave  even  before  it 
reaches  the  incident  shock  wave.  This  agrees  with  the  observation  made  by 
Bazhenova  and  Soloukhin  (1959)  that  the  merging  of  the  flame  front  with  the 
incident  shock  wave  is  not  a  necessary  condition  for  the  detonation  wave 
formation.  The  simulations  also  show  a  large  pressure  overshoot  when  the 
forward  moving  wave  coalesces  with  the  incident  shock  wave.  This  has  been 
observed  earlier  in  the  experiments  of  Urtiev  and  Oppenheim  (1966). 

The  numerical  simulations  presented  in  this  paper  show  that  the  one¬ 
dimensional  reactive  flow  model  with  detailed  chemical  kinetics  can  be  used 
to  elucidate  some  of  the  details  of  weak  ignition  behind  incident  shocks  and 
the  subsequent  transition  to  detonation.  The  predictions  of  the  model  are  in 
qualitative  agreement  with  experimental  observations.  However,  the  model 
does  not  include  multi-dimensional  phenomena  such  as  turbulence  and  boundary 
layer  growth  which  play  an  important  part  in  any  quantitative  study  of  the 
transiton  to  detonation.  Currently  a  two-dimensional  reactive  shock  model 
exists  but  it  uses  a  parameterized  model  for  energy  release  (Oran  et  al.  , 
1981b}  Oran  et  al. ,  1982).  The  model  is  now  being  extended  to  include  a 
detailed  chemical  kinetic  scheme.  Calculations  with  this  new  model  would 
show  the  effects  of  transverse  waves  and  boundary  layers  on  the  transition  to 


detonation: 


ACKII OVTLEDG  EMENTS 


The  authors  greatfully  acknowledge  suggestions 
conversations  with  Jay  Boris,  Ted  Young,  and  Antoni 
assistance  from  Mrs.  Fran  Rosenberg.  This  work  has 
Office  of  Raval  Research  through  the  Raval  Research 


from  and  useful 
Oppenheim  and  editorial 
been  supported  by  the 
laboratory . 


c* 


D 


C-16 


REFERENCES 


Bahn,  G.S.  {1968).  Reaction  Rate  Compilation  for  H-O-N  System,  Gordon  and 
Breach,  New  York. 

Baulch,  D.L. ,  Drysdale,  D.C.,  Horne,  D.C.,  and  Lloyd,  A.C.  (1972).  Evaluated 
Kinetic  data  for  High  Temperature  Reactions,  Vol.  1,  Butterworths, 
London. 

Bazhenova,  T.V. ,  and  Soloukhin,  R.I.  (1959).  Gas  Ignition  Behind  the  Shock 
Wave.  Seventh  Symposium  (International)  on  Combustion,  The  Combustion 
Institute,  Pittsburgh,  pp.  866-875* 

Boris,  J.P.  (1976).  Flux-Corrected  Transport  Modules  for  Solving  Generalized 
Continuity  Equations.  Naval  Research  Laboratory  Memorandum  Report  3237, 
Washington,  D.C. 

Boris,  J.P. ,  and  Book,  D.L.  (1976).  Solution  of  Continuity  Equations  by  the 
Method  of  Flux-Corrected  Transport.  In  Methods  of  Computational 
Physics,  Academic  Press,  New  York,  Vol.  16,  Chap.  11,  pp.  85-129. 

Burks,  T.L. ,  and  Oran,  E.S.  (1980).  A  Computational  Study  of  the  Chemical 
Kinetics  of  Hydrogen  Combustion.  Naval  Research  Laboratory  Memorandum 
Report  UU1»6,  Washington,  D.C. 

Cohen,  N. ,  and  Westberg,  K.R.'  (1979).  Data  Sheets.  The  Aerospace 
Corporation,  Los  Angeles. 

Edvards,  D.H.,  Thomas,  G.O. ,  and  Williams,  T.L.  (1981).  Initiation  of 

Detonation  by  Steady  Planar  Incident  Shock  waves.  Combust.  Fleme.  1*3, 
187. 

Hampson,  R.F. ,  and  Garvin,  D.  (1975)*  Chemical  ■Kinetic  and  Photochemical 

Data  for  Modeling  of  Atmospheric  Chemistry .  NBS  Technical  Note  866,  U.S. 
National  Bureau  of  Standards,  Washington,  D.C. 


Lee,  J.H.,  and  Ramamurthi,  K.  (1976).  On  the  Concept  of  the  Critical  Size  of 
a  Detonation  Kernal.  Combust.  Flame.  27,  331. 

Lee,  J.H.  (1977).  Initiation  of  Gaseous  Detonation.  Ann.  Pev.  Efrys.  Chem. 
28,  75- 

Lloyd,  A.C.  (197*0..  Evaluated  and  Estimated  Kinetic  Data  for  Gas  Fhase 

Reactions  of  the  Hydropero^l  Radical.  Int.  J.  Chem.  Kinetics.  6,  169. 

Meyer,  J.W. ,  and  Oppenheim,  A.K.  (1971a).  On  the  Shock-Induced  Ignition  of 
Explosive  Gases.  Thirteenth  Symposium  (International)  on  Combustion, 

The  Combustion  Institute,  Pittsburgh,  pp.  1153-1161*. 

Meyer,  J.W. ,  and  Oppenheim,  A.K.  (1971b).  Coherence  Theory  of  the  Strong 
Ignition  Limit.  Combust.  Flame.  17,  65* 

Olson,  D.B. ,  and  Gardiner,  W.C.  (1977).  An  Evaluation  of  Methane  Combustion 
Mechanisms.  J.  Phys.  Chem.  8l,  25lU. 

Oran,  E.S.  Young,  T.R. ,  and  Boris,  J.P.  (1979).  Application  of  Time- 
Dependent  Numerical  Methods  to  the  Description  of  Reactive  Shocks. 
Seventeenth  Symposium  (International)  on  Combustion.  The  Combustion 
Institute,  Pittsburgh,  pp.  1*3-5**. 

Oran,  E.S. ,  and  Boris,  J.P.  (1981a).  Detailed  Modelling  of  Combustion 
Systems.  Prog.  Energy  Combust.  Sol.  7,  1* 

Oran,  E.S. ,  and  Boris,  J.P.  (1981b).  Weak  and  Strong  Ignition-II.  Naval 
Research  Laboratory  Memorandum  Repot  1*671,  Washington,  D.C.  (also  to 
appear  in  Combust.  Flame). 

Oran,  E.S.,  Young,  T.R.,  Boris,  J.P. ,  and  Cohen,  A.  (1961a).  Weak  and  Strong 
Ignition-I.  Naval  Research  laboratory  Memorandum  Report  1*661*, 
Washington,  D.C.  (also  to  appear  in  Combust.  Flame) 


C- 18 


■y 


<• 


n 


v 


Oran,  E.S.,  Boris,  J.P.  ,  Young,  T. ,  Flanigan,  M. ,  Burks,  T. ,  and  Picone,  M. 
(I9?lfc).  Numerical  Siiailations  of  Detonations  in  Pydrogen-Air  and 
Methane-Air  Mixtures.  Eighteenth  Symposium  (International)  on 
Combustion,  The  Combustion  Institute,  Pittsburgh,  pp.  I6UI-16L9. 

Oran,  E.S.,  Young,  T.R.  Boris,  J.P. ,  and  Picone,  J.M.  (1962) .  A  Study  of 
Detonation  Structure:  The  Formation  of  threacted  Gas  Pockets. 

Presented  at  the  Nineteenth  Symposium  (International)  on  Combustion,  The 
Combustion  Institute,  Pittsburgh. 

Strehlow,  R.A. ,  Crooker,  A.J.,  and  Cusy,  R.E.  (1967).  Detonation  Initiation 
Behind  an  Accelerating  Shock  Wave.  Combust.  Flame.  11,  339. 

Stull,  D.R.,  and  Prophet,  H.  (Eds.)  (1971).  JANAF  Ihermochemlcal  tables. 

Rat.  Stand.  Ref.  Data  Serv. ,  Ho.  37,  2nd  edition,  U.S.  National  Bureau 
of  Standards,  Gaithersburg. 

Urtiew,  P. A. ,  and  Oppenheim,  A.K.  (1966).  Experimental  Observations  of  the 
Transition  to  Detonation  in  an  Explosive  Gas.  Proc.  Roy.  Soc.  A295,  13. 

Urtiew,  P.A. ,  and  Oppenheim,  A.K.  (1967).  Detonative  Ignition  Induced  by 
Shock  Merging.  Eleventh  Symposium  (International)  on  Combustion.  The 
Combustion  Institute,  Pittsburgh,  pp. 665-670. 

Williams,  F.A.  (1965).  Combustion  Theory.  Addison4Jesley,  Reading,  p.  2. • 

Young,  T.R. ,  and  Boris,  J.P.  (1977).  A  Numerical  Technique  for  Solving  Stiff 
Ordinary  Differential  Equations  Associated  with  the  Chemical  Kinetics  of 
Reactive-Flow  Problems.  J.  Phys.  Chem.  8l,  2k2k. 

Young,  T.R.  (1980).  CHEMEQ  -  A  Subroutine  for  Solving  Stiff  Ordinary 

Differential  Equations.  Raval  Research  Laboratory  Memorandum  Report 
1*091,  Washington,  D.C.  _ 


C- 19 


Ikble  I.  H2-02  Elementary  Reactive  Mechanism 


H  +  HO  X  0  +  H2 

l.l*0(-ll*) 

1.00 

3.50(+03) 

[1] 

3.00(-lM 

1.00 

l*.l*8(+03) 

[ll 

H  +  H02  ♦  H2  +  02 

l*.20(-ll) 

0.00 

3. 50(+02) 

[1] 

9.10(-11) 

0.00 

2.9l(+0l*) 

ID 

H  +  H02  t  HO  +  HO 

1*.20(-10) 

0.00 

9-50(+02) 

ll] 

2.00(-ll) 

0.00 

2.02(+0l*) 

111 

H  +  H02  Z  0  + 

8.30(-ll) 

0.00 

5. 00(+02) 

12] 

1.75(-12) 

0.1*5 

2.8U(+OU) 

kr  =  kf/K( 

H  +  H202  t  H02  +  H2 

2.80(-12) 

0.00 

1.90(+03) 

ll] 

1.20(-12) 

0.00 

9.1*0(+03) 

ll] 

H  +  H202  $  HO  +  H20 

5.28(-10) 

0.00 

l*.50(+03) 

ll] 

3.99(-10) 

0.00 

l*.05(+OU) 

kr  «  kf/K( 

HO  +  H2  $  H  +  H20 

1.83(-15) 

1.30 

1. 81*(  +03) 

[3] 

1.79(-lit) 

1.20 

9.6l(+03) 

[3] 

HO  +  HO  i  H2  +  02 

1.09(-13) 

0.26 

l.l*7(+Ol*) 

kf  *  krKc 

2.82(-ll) 

0.00 

2.1+2(1-01*) 

T\h]  C 

HO  +  HO  t  0  +  H20 

1.00(-l6) 

1.30 

0.00(+00) 

13] 

3.20(-15) 

1.16 

8.77(+03) 

kr  =  kf/K( 

HO  +  H02  ♦  H20  +  02 

8.30{-U) 

0.00 

5«03(+02) 

1 5 1 

2.38(-10) 

0.17 

3.69(+0l*) 

kr  =  VK« 

HO  +  H,0  $  HO,  +  H, 

1.70(-11) 

0.00 

9.10(+02) 

ll] 

l+.70(-ll) 

0.00 

1.65(+0l+) 

111 

HO  +  H2  ♦  HO  +  H20 

1.20(-12) 

0.00 

9.1*l(+03) 

ll*] 

1.33(-1U) 

0.1*3 

3.62(+0U) 

kr  -  VK« 

ho2  +  ho2  t  h2o2  +  o2~ 

3.00{-ll) 

0.00 

5.00(+02) 

12] 

1.57(-09) 

-0.38 

2.20(+0l*) 

kr  «  kf/K( 
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Reaction 

ki  =  AT5 

A(b) 

exp  (-C/T)(a) 

2 

c(b) 

References 

0+  HO  X  H  +  02 

2. 72(-12) 

0.28 

-8.10(+01) 

3.70(-10) 

0.00 

8.1*5(+03) 

f [i]  r  c 

0  +  H02  t  HO  +  02 

B.32(-ll) 

0.00 

5.03(+02) 

151 

2.20(-ll) 

0.18 

2.82(+0l*) 

kr  «  kf/Kc 

o  +  e2o2  i  h2o  +  o2 

1.  *»0( -12) 

0.00 

2.12(+03) 

12] 

5- 70(  — lU) 

0.52 
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1.1*0(  -12) 

0.00 

2.13(+03) 
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2.07(-15) 

0.61* 

8.23(+03) 

kr  =  kf/Kc 
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Ill 
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0.00 

l*.83(-Ol*) 
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6.20(-26) 

-2.00 
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HI 

5.80(-09) 

0.00 

5.29(+0l*) 

[1] 
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l».lU(-33) 
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HI 
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2.30(+0l*) 

[1] 

HO  +  HO  +  M  X  H202  ♦  M 

2.50(-33) 

0.00 

-2.55(+03) 

ID 

2.00(-07) 

0.00 

2.29(+0k) 

111 

O+H+MJHO+M 

8. 28( —29) 
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0.00(+00) 
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2.33(-10) 
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(a)  Bimolecular  reaction  rate  constants  are  given  in  units  of  cm3/ (molecule  sec). 

Te  molecular  reaction  rate  constants  are  given  in  units  of  cm6/ (molecule 2  sec). 

(b)  Exponentials  to  the  base  10  are  given  in  parenthesis;  i.e. ,  1.00(-10)  * 

1.00  *  10"l°. 

(c)  The  references  are:  (l)  Baulch  et  al. ,  1972;  (2)  Hampson  and  Garvin,  1975; 

(3)  Qohen  and  Westberg,  1979;  (**)  Olson  and  Gardiner,  1977;  (5)  Lloyd,  197**; 

(6)  Bahn,  1968. 


liable  II.  Parameters  for  the  Weak  Ignition  Study 


Fhrameter 


Driver  Section 


Length  30  cm 

Gas  Mixture  He 

Initial  Itenperature  298  K 

Initial  Pressure  9  atm 


Incident  Shock  Velocity 
Itenperature  behind  Incident  Shock 
Pressure  behind  Incident  Shock 
Chapman-Jouguet  Detonation  Velocity 


Driven  Section 


370  cm 

E2:02:N2/2:1:L 
298  K 
0.1  atm 
1.1*  xlO5  cm/s 
918  K 
1.39  atm 
1.9*105  cm/6 


Table  III.  Time  History  of  the  Forward  Moving  Reaction  Wave. 


Time  Velocity  of  Wave  Fluid  Velocity  Sound  Velocity 

(us)  (cm/s)  (cm/s)  (cm/s) 


Ahead 

Behind 

Ahead 

Behind 

609.03 

2.70(+05) 

1.09(+05) 

1.30(+05) 

7.15(+OU) 

1.10(+05) 

618.63 

2.72(+05) 

1.08(+05) 

1.50(+05) 

7.15(+0U) 

1.12(405) 

626.03 

2.75(+05) 

1.08(+05) 

1.60(+05) 

7.10(+0U) 

1.15(+05) 

631.83 

2.80(+05) 

1.08(+05) 

1.85(+05) 

7.08(+0U) 

l.lU(+05) 

663.83 

3.20(+05) 

1.07(+05) 

2.10(+05) 

7.20(+0U) 

l.l8(+05) 

688.13 

2.60(+05) 

0 

2.10(+05) 

U.0U(+0U) 

1.28(+05) 

Tkble  IV.  Tine  History  of  the  Backward  Moving  Reaction  Wave 


Time  Velocity  of  Wave  Fluid  Velocity  Sound  Velocity 

(us)  (cm/s)  (cm/s)  (cm/s) 


Ahead 

Behind 

Ahead 

Behind 

609.03 

•  -T.2(+0U) 

1.08(+05) 

8.3(+05) 

7.2(+0l0 

1.06(+05) 

618.63 

-7.3(+0U) 

1.08(+05) 

1.0(+0h) 

7.2(+0l*) 

l.lK+05) 

626.03 

-7. 5(+0l») 

1.09(+05) 

6.0(*0k) 

7.2(+0U) 

1.13(+05) 

631.83 

-7.6(+0l*) 

1.09(+05) 

5.5(+OU) 

7.2(+0U) 

l.lU(+05) 

FIGURE  CAPTIONS 


Figure  1.  A  schenatic  representation  of  the  adaptive  grid  used  for  the 
numerical  simulations. 

Figure  2.  A  position-time  diagram  of  the  main  events  occurring  in  the  shock 
tube  sisulation:  S,  incident  shock;  CS,  contact  surface;  RC, 
reactive  centers,  RW,  reaction  wave;  (a)  pre-ignition  regime,  (b) 
quasi-steady  shock -react  ion  complex,  (c)  formation  of  reactive 
centers,  (d)  hot  spot  formation  leading  to  overdriven  detonation 
and  (e)  detonation  relaxation. 

Figure  3.  Hie  spatial  variation  of  the  temperature  and  the  05  mole  fraction 
between  the  contact  surface  and  the  shock  at  135  Us. 

Figure  1*.  3he  spatial  distribution  of  the  tenperature  between  the  contact  • 
surface  and  the  shock  at  four  different  times. 

Figure  5*  The  spatial  variation  of  the  temperatue  and  the  OH  mole  fraction 

between  the  cntact  surface  and  the  shock  at  (a)  3^8  Us  and  (b)  1*1*5 
us. 

Figure  6.  Time  history  of  the  development  of  a  hot  spot. 

Figure  7.  A  position-time  diagram  shewing  the  ignition  of  flamelets  and  the 
transition  to  detonation:  A,  hot  spot  and  B,  flamelets. 

Figure  8.  The  pressure,  velocity  and  tenperature  distributions  across  the 
system  at  (a)  599*33  Us  and  (b)  609*03  us. 

Figure  9*  The  pressure  and  tenperature  distributions  across  the  system  at 
(a)  618.63  us,  (b)  626.03  Us  and  (c)  631.83  Us. 

Figure  10.  The  pressure  and  tenperature  distributions  across  the  system  at 
(a)  663.83  us  and  (b)  688.13  Us. 
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